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ABSTRACT 

In this paper, we propose a novel reduced-rank algorithm for direc- 
tion of arrival (DOA) estimation based on the minimum variance 
(MV) power spectral evaluation. It is suitable to DOA estimation 
with large arrays and can be applied to arbitrary array geometries. 
The proposed DOA estimation algorithm is formulated as a joint op- 
timization of a subspace projection matrix and an auxiliary reduced- 
rank parameter vector with respect to the MV and grid search. A 
constrained least squares method is employed to solve this joint op- 
timization problem for the output power over the grid. The pro- 
posed algorithm is described for problems of large number of users' 
direction finding with or without exact information of the number 
of sources, and does not require the singular value decomposition 
(SVD). The spatial smoothing (SS) technique is also employed in 
the proposed algorithm for dealing with correlated sources problem. 
Simulations are conducted with comparisons against existent algo- 
rithms to show the improved performance of the proposed algorithm 
in different scenarios. 

Index Terms — Direction of arrival (DOA) estimation, array 
processing, joint iterative methods, reduced-rank methods. 

1. INTRODUCTION 

Direction of arrival (DOA) estimation techniques have been widely 
employed in many fields related to array processing JT]. Numerous 
DOA estimation approaches have been considered to date. Among 
them are the Capon |2], the conventional subspace-based methods 
that require the singular value decomposition (SVD), such as MU- 
SIC (3) and ESPRIT Q, and more recent subspace techniques that 
do not require the SVD, such as the auxiliary vector (AV) estimation 
algorithm [5 1 and the conjugate gradient (CG) algorithm [6 |. 

The Capon DOA estimation method minimizes the output power 
of the undesired interferences while maintaining a constant gain 
along the look direction. By computing and plotting Capon's spec- 
trum over the possible scanning directions, the DOAs can be es- 
timated by locating the peaks in the spectrum. The estimation 
accuracy of the Capon method strongly depends on the number of 
snapshots and the array size. The subspace-based MUSIC and ES- 
PRIT algorithms exploit the eigen-structure of the input covariance 
matrix to decompose the observation space into a signal subspace 
and a corresponding orthogonal noise subspace. ESPRIT has bet- 
ter performance by employing a displacement invariance in some 
specific array structures. The developed eigen-decomposition algo- 
rithms are described in Q, (§]. The previously reported methods 
suffer from correlated sources and high computational complexity 
due to the eigen-decomposition procedure. The AV and CG esti- 



mation algorithms were proposed recently. The AV method is de- 
veloped based on the orthogonality of an extended non-eigenvector 
signal subspace with the true signal subspace and the scanning vec- 
tor itself. As the scanning vector drops in the signal subspace, the 
DOAs are determined by finding the collapse in the extended signal 
subspace. The CG method can be considered as an extended version 
of the AV method since it applies the residual vectors in place of the 
AV basis. Both algorithms show dominate in severe conditions with 
a small number of snapshots and at low SNR for both correlated and 
uncorrelated sources. However, they work inefficiently with a large 
number of sources or without exact information about the number 
of sources beforehand. 

In this paper, we propose a DOA estimation algorithm by em- 
ploying a novel reduced-rank signal processing strategy. The pro- 
posed algorithm is based on a joint iterative subspace optimization 
(JISO) and grid search with respect to the MV power spectrum eval- 
uation. The implementation of the proposed DOA estimation al- 
gorithm amounts to designing a subspace projection matrix and an 
auxiliary reduced-rank parameter vector with respect to the MV cri- 
terion. We present a constrained least squares algorithm for jointly 
estimating the subspace projection matrix and the auxiliary reduced- 
rank parameter vector that calculate the output power over the pos- 
sible scanning directions. The proposed algorithm is more practical, 
in comparison with the existing algorithms, since it is not limited by 
the array structure, does not require the SVD procedure, and works 
without information of the source number, which will be shown in 
simulations. The estimation accuracy is also satisfied with a large 
number of sources' direction finding. In addition, the spatial smooth- 
ing (SS) technique, which was devised by Evans in |9| and further 
developed by Shan in 1101 , is employed in the proposed DOA esti- 
mation algorithm for dealing with the problem caused by correlated 
sources. 

The rest of this paper is organized as follows: we outline a sys- 
tem model for DOA estimation and present the problem statement in 
Section 2. Section 3 derives the proposed DOA estimation algorithm 
and analyzes the complexity. The application of the SS technique in 
the proposed algorithm is also introduced briefly in this part. Simula- 
tion results are provided and discussed in Section 4, and conclusions 
are drawn in Section 5. 

2. SYSTEM MODEL AND PROBLEM STATEMENT 
2.1. System Model 

Let us suppose that q narrowband signals impinge on a uniform lin- 
ear array (ULA) of m (m > q) sensor elements. Note that the pro- 
posed DOA estimation algorithm can be applied to arbitrary array 



structures. An extension to arbitrary arrays will be sought in a fu- 
ture work. The ULA here is adopted for using the SS technique and 
reaching a fair comparison with ESPRIT, which is applied to some 
specific array structures. The ith snapshot's vector of sensor array 
outputs x(i) G Q mxl can be modeled as 



x(i) — A(0)s(i) + n(i), i 



(1) 



where = [9 , 6> 9 -i] T G C qxl is the signal DOAs, A(0) = 
[a(6o), ■ a(9 q -i)] G C mxq is the matrix that contains the signal 
direction vectors a(9k), where a(9k) = [1, e~ 2n -'^ cosek , . . . , 
e -2„( m -i)£ cosefc] T g c mxl , (k = 0, . . . , q - 1), A c is the 
wavelength, and d (d — A c /2 in general) is the inter-element dis- 
tance of the ULA. To avoid mathematical ambiguities, the direc- 
tion vectors a(9k) are considered to be linearly independent (§]. 
s(i) G 1Z qxl is the source data. n(i) G C mxl is the white sensor 
noise, which is assumed to be a zero-mean spatially and Gaussian 
process, TV is the number of snapshots, and (-) T denotes transpose. 

2.2. Problem statement 

Based on the MV output power spectrum (or Capon output power 
spectrum), [2J, 1111 . the output power to each scanning direction for 
DOA estimation is expressed by 



= argminiUfl Rwg 



subject to w e a(9) 



(2) 



where 9 is the estimated direction and wg — [we,i, ■ ■ ■ , U!g ym ] T G 
C mxl is the weight vector corresponding to the current scanning 
direction 9. (-) H denotes Hermitian transpose. R is the data covari- 
ance matrix 



R = E{x(i)x H (i)] = A(0)R s A H (0) + , 



(3) 



where R s = E[s(i)s H (i)] denotes the signal covariance ma- 
trix, which is diagonal if the sources are uncorrelated and is 
nondiagonal and nonsingular for partially correlated sources, and 



E[n(i)n l 



= a n I with I being the corresponding identity 



matrix. 

The MV (Capon) power spectrum estimation algorithm attempts 
to minimize the contribution of the total output power while main- 
taining an unity gain along a look direction 9. By optimizing the 
weight vector wg and obtaining the output power for all possible di- 
rections 9 G (0°, 180°), the DOAs can be determined by finding the 
peaks in the output power spectrum. The weight solution is (2), 1111 



Wg 



R- 1 a(9) 
a H (e)R- 1 a(9) 



(4) 



Substituting Q into l|2}, DOA estimation based on the MV 
(Capon) power spectrum is given by 



9mv = arg max [a H (9)R 1 a(9)] 1 



(5) 



Note that complete knowledge of R cannot be obtained in prac- 
tice. We may use a sample-average recursion to estimate this input 
covariance matrix, which is given by 



1 N 



(6) 



Where R is not invertible if the number of available snapshots is 
less than the number of sensors (TV < m). It can be implemented by 
employing the diagonal loading technique QD- 

The above MV based DOA estimation method suffers from a 
heavy computational load for large m due to the matrix inversion 
and works inefficiently in the presence of correlated sources. Fur- 
thermore, the performance is inferior when large number of sources 
appear in the system. 

3. PROPOSED DOA ESTIMATION ALGORITHM 

In this section, we employ a reduced-rank strategy to perform DOA 
estimation. This is carried out via the proposed joint iterative sub- 
space optimization (JISO) according to the MV criterion for estimat- 
ing the subspace projection matrix and the auxiliary reduced-rank 
parameter vector followed by a grid search. 

3.1. Proposed Reduced-Rank DOA Estimation Scheme 

We introduce a subspace projection matrix TV = [ti, ta, • • • , t r ] 
G C mxr , which is responsible for the dimensionality reduction, to 
project the mxl received vector x(i) onto a lower dimension, yield- 
ing 

x(i) = T"x(i) (7) 

where t, = [ti,i,t 2 , t , ■ ■ ■ ,tm/ G C mxl , (I = 1, . . . ,r) makes 
up the subspace projection matrix TV, x(i) G C rxl is the projected 
received vector, and in what follows, all r dimensional quantities are 
denoted with a "bar", r < m is the rank. An auxiliary filter with the 
reduced-rank weight vector f g = [fe,i, fe,2, ■ ■ ■ , fe,r] T G C xl 
is applied after the projection procedure. The aim of TV is to ex- 
tract the key features of the original input vector x(i) and form the 
reduced-rank input vector x(i). The auxiliary reduced-rank weight 
vector fg works on x(i) for obtaining the output power with respect 
to the current scanning direction 9. Since the procedure is operated 
with a lower dimension r, the computational complexity will be re- 
duced if r « m. Since DOA estimation depends on the number of 
sensor elements m and on the eigenvalue spread of the input covari- 
ance matrix, the proposed reduced-rank estimation scheme will ex- 
hibit improved performance under conditions where m is large 11121 . 
Following the MV DOA estimation in {2j, the proposed optimization 
problem can be expressed by 



9jiso = axg min f g T^RT r f e 



subject to f H T"a{9) = 1 



(8) 



We find that the minimization with respect to ([8} is equivalent to 
the joint optimization of the subspace projection matrix T r and the 
auxiliary reduced-rank weight vector f g . After obtaining T r and 
fg, DOA estimation can be determined by plotting the output power 
spectrum for the possible directions and searching for peaks that cor- 
respond to the DOAs of the sources. It is worth noting that, for 
r — 1, the novel scheme becomes a conventional full-rank MV 
scheme with an additional weight parameter fg that provides an am- 
plitude gain. For r > 1, the signal processing tasks are changed and 
T r and f g are optimized for obtaining the proposed output power 
spectrum for the possible directions. 

3.2. Proposed Joint Iterative Subspace Optimization Algorithm 

The challenge left to us is how to efficiently compute the subspace 
projection matrix T r and the auxiliary reduced-rank weight vector 



fg for solving the optimization problem ((8}. We propose a con- 
strained least squares (LS) algorithm to solve this joint optimization 
problem. The constraint in (|8j can be incorporated by the method of 
Lagrange multipliers 11131 in the form 

•7 = i2 ai ~ l \fe H {^(i)x(l)\ 2 + X[f e H (i)T^[t)a(e) - l] 

(9) 

where a is a forgetting factor, which is a positive constant close to, 
but less than 1, and A is a scalar Lagrange multiplier. Fixing fg(i), 
computing the gradient of ((9} with respect to T r (i), yields 

i 



= R(i)T r (i)f 9 (i)f 9 H (i) + \ Tr a(0)f g H (i) 



(10) 



where R{i) = J^Li a?- l x(l)x H (I) € C mxm is the estimated co- 
variance matrix According to [13], R(i) can be written in a recursive 
form as 

R(i) = afl(t - 1) + x(i)x"(i) (11) 

Making the above gradient terms equal to zero, multiplying 
f g (z) from the right of both sides, and rearranging the expression, it 
becomes, 

T r {i)f g = -\ Tr R\i)a{e) (12) 

where R (i) is invertible by employing the diagonal loading tech- 
nique. 

If we define p(i) = —\T r R (i)a(6), the solution of T r (i) 
can be regarded to find the solution to the linear equation 



T r (i)f e =P(i) 



(13) 



In order to find an unique solution for T r (i), we express the 
quantities involved in J 1 3b by 



T r (i) = 



(14) 

The problem in d!3t is equivalent to find p^ (i) (j = 1, . . . , m) 
for satisfying 



"PxW 




7«,i(i)" 




'pi(i)~ 


p 2 (*) 








P2(i) 




; fS) = 












JeA*). 







min Wpj&W , subject to p d (i)f e =pj(i) 
which is obtained by using the Lagrange multiplier method 

fe H (i) 



ll/ fl (0ll s 



and thus the projection matrix is 



T r {i) = P (i) 



Substituting the definition of p(i) into J 17t . we have 



T r {i) = -\ Tr R (i)a(6)- 



\f e (i) 



(15) 



(16) 



(17) 



(18) 



Table 1. The JISO algorithm for each scanning direction 
Initialization: 

TY(0) = [iT o^ x(m _ r) ] 

f e (0) = (T?(0)a(8 n ))/(\\T?(0)a(e n )\\ 2 ) 
Update for each time instant i = 1, . . . , N 

x(i) = Tf (i - l)x(i) 

a(6 n ) = T?(i-l)a(6 n ) 

R(i) = aR(i - 1) + x(i)x H (i) 

R(i) = aR(i - 1) + x(i)x H (i) 

f g (i) = R- 1 (i)a(e)/(a H (0)R- 1 (i)a(e)) 

T (j \ _ R 1 (i)o.(0) fe H (i) 

Output power 

PnsojOn) - l/(a g (e n )R- 1 a(e n )) 



The multiplier \r r can be solved by incorporating ( 112) with the 
constraint in ([§}, which is 



a ff (0)H \i)a(6) 
Substituting dl9| l into d!8t , we get the projection matrix 

R-\i)a(6) f e H (i) 



(19) 



T r (i) 



a H {0)R \i)a{9) 11/eW 



(20) 



At the same time, fixing T r (i), taking the gradient of (|9j with 
respect to fg(i), and making it equal to a null vector, we obtain 



= fl(*)/ fl (») + A A T?(*)a(0) 



(21) 



where = E)!=i a ' l x(l)x H (l) G C rxr is the estimate 

of the reduced-rank covariance matrix R = E[x(i)x H (i)] — 
T?E[x(i)x H (i)]T r . 

Following the same procedures for calculating T r (i), we obtain 
the result for the auxiliary reduced-rank weight vector f g (i) 



fe® = -\ fe R-\i)T?{i)a{0) 



d H (d)R- 1 (i)a(e) 
H- 1 (i)o(^) 



(22) 
(23) 

(24) 



a H (0)R- 1 (i)a(6) 

where a(6) — T^a(6) g C x 1 is the projected steering vector with 
respect to the current scanning direction. Note that (1241 is similar in 
form to © if we do not consider the time instant i. The proposed 
reduced-rank weight vector f g (i) is more general when dealing with 
DOA estimation, namely, for r = m, it is equivalent to the MV 
weight vector, and, for 1 < r < m, it operates under lower dimen- 
sions for reducing the complexity and improving the performance. 



3.3. DOA Estimation 

After snapshots, substituting the weight solution f g expressed in 
J241 with respect to the possible scanning directions 9 £ (0°, 180°), 
and the subspace projection matrix TV in i20l into JgJ, we obtain the 
corresponding output power spectrum for DOA estimation 

Pjiso(0„) = (a H (^ n )l?- 1 a(e n ))~ 1 (25) 

where the scanning direction 6 n — nA°, A° is the search step, and 
n = 1,2,..., 180°/A°. For a simple and convenient search, we 
make 180°/A° an integer. We use a similar form to that of (II 11 

for estimating R. The proposed JISO algorithm for each scanning 
direction 6 n is summarized in Table [TJ where T r (0) and f g (0) are 
initialized to ensure the constraint. The proposed algorithm provides 
an iterative exchange of information between the projection matrix 
and the reduced-rank weight vector, which leads to the improved 
performance. 

The output power in ((25} is much higher if the scanning direction 
n — Ok, (k = 0, . . . , q — 1), which corresponds to the transmit- 
ted sources, compared with other scanning angles that correspond to 
the noise level. Therefore, the output power spectrum shows peaks 
with respect to the sources when we plot it through the whole search 
range. 

Considering correlated sources, we can use the SS technique |9| 
in our proposed algorithm. It is based on averaging the covariance 
matrix of identical overlapping arrays and so requires an array of 
identical elements equipped with some form of periodic structure. 
We divide the ULA into overlapping subarrays of size n, with ele- 
ments {1, . . . , n} forming the first subarray, elements {2, . . . , n+1} 
forming the second subarray, etc., and J = ra — n+1 as the number 
of subarrays. Note that the selection of n needs to follow n > q and 
J > q 1 10 1 . The SS preprocessing scheme operates on the input vec- 
tor x(i) to obtain each subarray vector Xj(i), where j — 1, . . . , J. 
The proposed JISO algorithm is followed for DOA estimation. We 
denominate this SS-based algorithm as JISO-SS, which is summa- 
rized in TableU where TV, SS € C nxr and f e ss £ C rxl are the SS- 
based subspace projection matrix and auxiliary reduced-rank weight 
vector, respectively. A ss = [a ss (8 ), ■ ■ ■ , a ss (8 q -i)] £ C nxq is the 
matrix that contains the direction vectors a ss (9k), where ^(Ok) = 

[l,e- 2 ^^cos0 fe ,..., e - 2 ^ ( "- 1) *cose fc ] T £ C nxl ,(k = 

0, . . . , q - 1), D = diagje- 2 ^ £ cos e ° , . . . , e~ 27rj '£ cos e "- 1 } 
£ C qxq , and the subscript "ss" denotes that it is for the SS-based pro- 
posed algorithm. aSj lSS (i) £ C n x 1 is the input vector at the jth subar- 
ray, Xj,ss(i) £ C rxl is the corresponding reduced-rank input vector, 
and n JiSS (i) £ C nxl is the white sensor noise. Pj, ss (i) £ C nxn 
and Pj ttis (i) £ C rxr are the full-rank and reduced-rank covariance 
matrices of the jth subarray, respectively, at time instant i. Since it 
is a well-known technique for dealing with correlated sources, the 
details have been omitted but related references can be found in (9), 

ED. 

Checking the computational complexity, the conventional 
Capon 1 2 1, MUSIC (3) and ESPRIT algorithms work with 
0(m 3 ), and the recent AV (5) and CG algorithms have a higher 
computational cost [6]. With respect to the proposed algorithm, 
R (i) costs 0(m 3 ) but is invariable for the grid search, namely, 
the result obtained for the first scanning direction can be used for the 
rest. The complexity of the proposed JISO algorithm for each itera- 
tion is 0(r 3 ), which is less complex than the AV or CG methods if 
r « rn for large arrays. The complexity of the proposed JISO al- 
gorithm with constrained LS optimization method is slightly higher 



Table 2. The JISO-SS algorithm for each scanning direction 
Initialization: 

TV lSS (0) = [I r O rx (n-r)] 

/ 6 , a (0) = (T« s (0)a ss (^))/(||T« s (0)a ss (e n )|| 2 ) 
Update for each time instant i = 1, . . . , N 
for j = 1, . . . , J 

asj, S8 (i) = A ss Z) J_1 s(i) + n,j, a (i) 
x J> (i)=T* s (i-l)a 5 j, 8S (i) 
Pj, ss (i) = x jiSS (i)xf^(i) 

^,-(0 = *i,-(0*A»(0 

end 

a ss (0„) = T^ ss (i - l)a 88 (0„) 

= oRj; - i) + i >;•; , p jim (i) 

i(i)=ai(i-i) + iE; =1 p i)S ,(i) 
f g>ss (i) = ic 1 (i)fi«(0»)/(a?(0n)£; 1 (»)s«(0n)) 

' ss( ' <(e„)H s ; 1 (i)a ss (9„) ilT^WF 
Output power 

PnsoMOu) = l/{aS (flrQfls^Mfln)) 



than the MUSIC and ESPRIT methods and lower than the AV and 
CG algorithms. The JISO-SS algorithm is marginally more complex 
than the JISO one due to the SS preprocessing. Actually, we can 
employ other methods to solve the joint optimization problem (e.g., 
stochastic gradient, recursive least squares B131 ) to avoid the matrix 
inversion for complexity reduction, which will be analyzed in the 
near future. 

4. SIMULATIONS 

Simulations are performed for an ULA with half wavelength in- 
terelement spacing. We compare the proposed algorithm with the 
Capon, MUSIC, ESPRIT, AV, and CG methods, and run K = 1000 
iterations to get each curve. The SS technique is employed for each 
algorithm to improve the performance. In all experiments, the BPSK 
signals' power is a'j — 1 and the noise is spatially and temporally 
white Gaussian. The search step is A° = 1°. The DOAs are consid- 
ered to be resolved if |#jiso — 0k\ < 

In Fig. Q] we consider the presence of q = 2 highly correlated 
sources separated by 3° with correlation value c = 0.9, which are 
generated as follows: 

si ~ A/"(0, <t 2 ) and S2 — csi + \J\ — c 2 S3 (26) 

where S3 ~ Af(0, a'i). The sensor elements number is m = 30 and 
input SNR = -2dB. We set the forgetting factor a = 0.998, the 
reduced dimension r = 6, and the diagonal loading 8 = 5 x 10~ 4 
for the covariance matrix inverse in (1181 and ( 1191 . The probability 
of resolution (5), HJ is plotted against the number of snapshots. The 
proposed algorithm outperforms other existing methods with small 
number of snapshots. The curves between the proposed and the MU- 
SIC algorithms are shown to intersect when the number of snapshots 
increases. The performance of the AV and CG methods can be seen 
to be inferior when compared to the proposed algorithm for all ob- 
servation periods. Regarding the SS-based algorithms, we set the 
subarray size ton = 26, which accords with [ 10] and reaches a high 
probability of resolution. The performance of the algorithms with 
the SS technique is improved and the proposed algorithm still has 
better performance than the existing ones. 




Number of Snapshots 



In the last experiment, we assess the performance of the pro- 
posed and analyzed algorithms with an uncorrect number of sources 
qw 7^ q known by the receiver. This is more practical since the exact 
sources number has to be determined by procedures with extra com- 
putation cost and time. We keep the scenario as that in Fig. [2] but 
assume an uncorrect number of sources q w — 9 instead of q — 10, 
and increase the number of snapshots for the operation. The fixed in- 
put SNR = OdB. In Fig. [3] the MUSIC and its SS-based algorithms 
start to work with large number of snapshots, and the ESPRIT and 
its SS-based algorithms fail to resolve DOA estimation with the in- 
crease of the snapshots since q is critical to the eigendecomposition 
for the partition of the signal subspace and the noise subspace in the 
input covariance matrix. Also, the design of the AV basis and CG 
residual vectors depends strongly on q. The Capon and its SS-based 
algorithms work well under this condition since they are insensitive 
to the number of sources. The same holds for the proposed and its 
SS-based algorithms, but both exhibit better performance and lower 
complexity. We consider q u , > q condition and get the same result. 



Fig. 1. Probability of resolution versus number of snapshots (sep- 
aration 3°, SNR= -2dB, q = 2, c= 0.9, m = 30, r = 6, 
S = 5 x 10~ 4 , a = 0.998, n = 26) 



Next, we consider the sources to be uncorrelated but increase 
the number of sources by setting q — 10. The input SNR = — 5dB 
and the number of sensor elements is set to m — 50. As can be 
seen in Fig. [2] the AV and CG methods are unable to obtain a DOA 
estimate with a large number of sources. The proposed algorithm 
demonstrates an improved performance and is the first to reach the 
highest resolution, as compared with the conventional Capon and 
the subspace-based MUSIC and ESPRIT methods, following the in- 
crease of number of snapshots. The subarray size is n = 41 in this 
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Fig. 2. Probability of resolution versus number of snapshots (sepa- 
ration 3°, SNR= -5dB, q= 10, m = 50, r = 6, 5 = 5 x 10~ 4 , 
a = 0.998, n = 41). 
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Fig. 3. Probability of resolution versus number of snapshots (sepa- 
ration 3°, SNR= OdB, q w = 9, m = 50, r = 6, 5 = 5 x 10~ 4 , 
a = 0.998, n = 41). 



5. CONCLUDING REMARKS 

We proposed a novel reduced-rank strategy to implement joint itera- 
tive subspace optimization and grid search for DOA estimation. The 
DOA estimation problem is formulated as a reduced-rank MV opti- 
mization problem. A subspace projection matrix is introduced to ob- 
tain the covariance matrix processed in the lower dimension so that 
computation cost is reduced and performance improved. An auxil- 
iary reduced-rank parameter vector is combined to realize the joint 
iterative optimization with respect to the MV output power for each 
scanning direction. By searching the possible directions, the DOAs 
can be determined by finding the peaks in the output power spec- 
trum. The proposed DOA estimation algorithm demonstrates ad- 
vantages under large array condition with uncorrelated or correlated 
sources. Its performance is not significantly influenced by some pa- 
rameters (e.g., the number of sources). In future work, we will pro- 



vide the analysis including the Cramer-Rao Bound (CRB) and com- 
pare it with the estimation accuracy of the proposed algorithm. We 
will also consider unitary versions of the proposed reduced-rank al- 
gorithm for ULA that do not reuqire grid search. 
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